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ABSTRACT 

The first non-trivial cumulant correlator of the galaxy density field Q21 is exam- 
ined from the point of view of biasing. It is shown that to leading order it depends on 
two biasing parameters b, and 62, and on 521, the underlying cumulant correlator of 
the mass. As the skewness Q3 has analogous properties, the slope of the correlation 
function —7, Q 3 , and Q 2 i uniquely determine the bias parameter on a particular scale 
to be b = 7/6(Q2i — Qs), when working in the context of gravitational instability 
with Gaussian initial conditions. Thus on large scales, easily accessible with the fu- 
ture Sloan Digital Sky Survey and the 2 Degree Field Survey, it will be possible to 
extract 6, and 62 from simple counts in cells measurements. Moreover, the higher or- 
der cumulants, Qn, successively determine the higher order biasing parameters. From 
these it is possible to predict higher order cumulant correlators as well. Comparing 
the predictions with the measurements will provide internal consistency checks on the 
validity of the assumptions in the theory, most notably perturbation theory of the 
growth of fluctuations by gravity and Gaussian initial conditions. Since the method 
is insensitive to f2, it can be successfully combined with results from velocity fields, 
which determine Q 6 /b, to measure the total density parameter in the universe. 
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1 INTRODUCTION 

Galaxies are likely to be biased tracers of the underlying 
density field, since surveys with different selection properties 
are con elated differently, i.e. they are biased with respect to 
each other. The amplitude ot the two point correlation tunc- 
tion is ambiguous and cannot be compared directly with 
measurements of the cosmic microwave background fluctua- 
tions, or dark matter simulations. Although the galaxy den- 
sity field can be a complicated function of the underlying 
mass field, on large scales the fluctuations are small and 
only the leading order contributions are important. Then 
the two -point function of galaxies becomes b 2 times the un- 



function, and §3, and Q3 are the skewness of the mass and 
galaxy field, respectively. The slope of the two-point func- 
tion determines the theoret i cal 173 using perturbation the 
ory (Peebles 198C[; |Fry 1994 ; Bernardcau li)92| : Juszkiewicz 



Bouchct, & Colombi 1993; Bernardcau 1994^ ~Unfortunately 



Q3 alone cannot yield both bias parameters, therefore ad- 
ditional information is needed. Note that adding Q4 does 
not help, since to leading order it depends on an additional 
bias parameter, 63, and at higher orders the situation is ex- 
actly analogous. Several attempts have been made in the 
literature to supply the missing piece of information an d to 
const rain the bias parameter by using the bi-spectrum ( ^ry 



derlyinj ; two-point correlation function, where b = b\ is the 
first coefficient in the Taylor expansion of the generaf (fo- 
cal) bias function, the galaxy density field as function of the 
underlying mass field. 

Higher order statistics can be used to alleviate this am- 
biguity by constraining the unknown value of b. The cumu- 
lants of the galaxy field are related to the cumulants of the 
underlying mass field and the biasing parameters (Fry and 
Gaztahaga 1993, hereafter FG93). For instance the skewness 
of the galaxy field was shown to be Q3 = 93/6 + 62/b 2 , where 
62 is the second order coefficient in the expansion of the bias 



1994 ; Matarrese, Verde, & Heavens 1997 



Scoccimarro et al. 



199£ ), configuration dependence of the three-point function 



( Fricman fc Gaztanaga 199q), and the skewness in elongated 
cells (Scoccimarro, Szapudi, & Frieman 1995) This Letter 
proposes a simple, robust alternative using the fir st non- 
trivial cumulant correlator ( Szapudi fc Szalay 1997a ) to pro- 
vide the needed additional constraint. This technique yields 
b via counts in cells measurements, thus it is possibly one 
of the simplest and most unique methods ever put forward. 
Note that there are alternatives besides higher order statis- 
tics for the determina tion of the bias, mo st importantly the 
use of velocity fields (Willick et al. 1997), which, however, 
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determine the combination of f2 0,6 /&. The simple method 
proposed here is insensitive to the cosmological parameter 
Q. Thus in combination with results from velocity flows the 
total density in the universe can be determined. 

Section §2 introduces the theory of biasing for cumu- 
lant correlators, §3 details how it can be applied to measure 
the bias parameters §4 discusses practicality, possible exten- 
sions, and limitations of the results. 



2 BIASING OF CUMULANT CORRELATORS 

Let us consider the moments and joint moments of the 
(smoothed) fluctuation field S, and denote the connected 
moments of the field with () c . The cumulants, gjv, can 
be defined by the moments as £ = (8 2 ), and (8 } c = 
qjviV^ -2 ^ -1 (N > 3). Similarly the cumulant correlators 
are defined through the joint moments = (<5i<52) and 
(6? 6$' ) c = qNMN N - l M M - 1 £ M+N - 2 £i (N + M > 3). The 
former is a well established tool for characterization of the 



non-Gaussianity of the galaxy distribution (Peebles 1980) 



while the latter were recently introduced by ( Szapudi fc Sza- 



Ia^L997a ) based on factorial moment correlators ( fennul 



3zalav_&_13oschanl 992 
3zapu*dr*eT^I!*T995^ ' 



Meiksin, Szapudi, & Szalay 1992 



If the galaxy density field is a general function of the 
underlying density field 8 g = f(S), similar quantities can 
be defined for the galaxies. Note that this Letter uses con- 
tinuum limit throughout, even though galaxies constitute a 
discrete process. This discreteness can be simply taken into 



accoun t by the use of factorial moments (Szapudi fe Szalay 



1993a) under the a ssumption of infinitesimal Poisson sam- 



pling ( Peebles 1980 ), and will not be mentioned further. On 
large scales, where the mass fluctuations are small, it makes 
sense to Taylor expand the galaxy density field in terms of 
the mass field (bias function) /(<$) = Y2k>o bk8 k /k\, where 
bk are called the fcth order bias parameters. FG93 have 
shown that to leading order it is possible to calculate the 
cumulants of the galaxy field in terms of cumulants of the 
underlying density field, and the bias parameters. Next their 
results are recapitulated briefly. For N < 3 the connected 
moments of the galaxy field are 



£ 9 



b 2 t 

qt/b + bi/b 2 , 



(1) 



where 6i = 6, the usual bias parameter for simplicity. The 
bias can be systematically calculated from the generating 
function log(e^ s ^ x ) , and explicit results up to sixth order 
can be found in FG93. The parameter bo is irrelevant for 
the higher order cumulants, it can be fixed by requiring 
(f(8)) = 0. The hierarchy is approximately inherited by the 
biased field to lowest non-vanishing order, even for multi- 
point estimators. Note that several other general calcula- 
tions were done in t he past for the correlations of biased 
fields, among others ( Kaiser 1984 ; Bardccn et al. 1986[ Grin 
3tein, & Wise 198(: ; |Ma tarrese, Lucchi n, fc Bonometto 1986^ 



Bzalay 1988| ; pzapudi 1994| ; |Matsubara 199 



Analogously to FG93, biasing of the cumulant cor- 
relators can be ca lculated from the gener ating function 
log ( e /(*i)*+/(*2)^ ( |Szapudi fc Szalay 1993a| ). It is a triv- 
ial although tedious matter to derive explicit formulae from 
it. Here only those results are shown which are needed for 



the applications in the next section, 
Q21 = q 2 i/b + b 2 /b 2 , 



(2) 



where the notation is analogous to the previous equation. 
Note that the dependence on the bias parameters is ex- 
tremely simple and exactly analogous to the corresponding 
results for Q3. This facilitates the solution of the equations 
in the next section. The higher order results will be shown 
explicitly elsewhere. Note that the generating function can 
be trivially generalized to higher than two-point estimators 
as well. 



3 MEASURING THE BIAS PARAMETER 

For models of structure formation where Gaussian initial 
conditions were amplified via gravity, the above formulae can 
be combined with perturbation theory results in the weakly 
non- linear regime ( |Peebles 1980| ; |Fry 1984fc Bernardcau 1994| ; 
Bernardeau 1996 ) 

q 3 = 34/21 - 7/3 (3) 
q 2 i = 34/21 - 7/6, 



where —7 is the slope of the mass correlation function at 
the particular scale considered. Note that if the bias is a 
slowly varying function of scale, the slope of the mass and 
galaxy correlation functions are the same on a given scale. 
This plausible condition is assumed for the rest of the pa- 
per. The above equations were found to be goo d approxima- 
tions in both galaxy c atalogs and simu l ations ( Peebles 198( ; 



Coles fc Frcnk 199l|; iGaztafiaga 1992 



Gaztafiaga 1994 ; Colombi et al. 1995 



Bouchet et al. 199£ 



Szapudi, Meiksin, & 



Nichol 1996| ; |Szapudi et al. 1998a| ). 

Thus given 7, Q3, and Q21, the equations can be solved 
for the bias parameters b, 62 • The result is 



7 



6(Q 2 







bo = 



7<32i(147 - 68) + 7Q 3 (68 - 7 7 ) 



a 



(4) 



(5) 



252(Q 2i 

These equations constitute the main result of this Letter. 
They become degenerate when 521 = f?3, i.e. 7 = 0. For typ- 
ical CDM spectra and for b 7^ the solution is well behaved 
on large scales, where the assumptions of this calculation 
hold. The above solution provides a simple and unique way 
of determining the bias parameter b. Note that Q21 and the 
Qn's up to a certain order fully resolve the degeneracy, as 
the Qjv's depend on N — 1 bias parameters, and 7, thus the 
6jv_i's can be successively calculated. From these further 
cumulant correlators can be predicted and compared with 
measurements, providing a set of consistency checks for the 
theory. 



4 DISCUSSION 

The previous section showed, that if —7, the slope of the 
correlation function, Q3, the third order cumulant, and Q21 
the first non-trivial cumulant correlator are measured, the 
bias parameter b can be calculated. Next the practicality of 
such a measurement in the near future will be reviewed. 
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According to (Colombi, Szapudi, & Szalay 1998), the 



third and fourth moments of the galaxy distribution will be 
measured with less than a few %, and 10% error, respec- 
tively, in the Sloan Digital Sky Survey (SDSS) on scales 
between 1 — 50/i -1 Mpc. Errors of the the same order can 
be expected for the future 2 Degree Field Survey (2dF), 
perhaps in a slightly narrower dynamic range of scales be- 
cause of its smaller size. Similar accuracy can be achieved 
for Q3 and Q21 , and even better accuracy for the correlation 
function on these scales. The previous considerations require 
that S « 1 both because of the Taylor expansion of the 
bias function, and reliance on perturbation theory. This still 
leaves quite a large dynamic range, roughly 10 — 50/i _1 Mpc, 
in which the bias parameters b, and 63, and their scale de- 
pendence can be determined accurately. Similarly, at higher 
orders using Qn, 6jv-i can be calculated successively. Thus 
on large scales the bias will be accurately measured facili- 
tating the comparison with the fluctuations in the cosmic 
microwave background and iV-body simulations. Note that 
the assumption of large scales is important for both the ap- 
plicability of the perturbation theory and Taylor expansion. 
On small scales, where the hierarchical assumption applies, 
i.e. the t hree-point fu nction ( = (53(^1^2 + £2^3 + £3£i) f° r 
N = 3 (|Peebles 198C|), Q21 - Q3 = &/2, as found in the 
APM (Szapudi & Szalay 1997a). This leads to an asymptot- 



doubt that this simple method will give a unique determi- 
nation of the bias on large scales from the SDSS, and from 
the 2dF. 

There are several further developments possible for the 
simple theory of the previous sections. First of all, it is en- 
tirely trivial to calculate the effects of bias on higher or- 
der cumulant correlators, given the generating function of 
§1. Similarly, the generating function generalizes easily for 
multipoint statistics. Projected catalogs can provide another 
way of estimating the needed three quantities, although de- 
projecti on in the weakly non-linear regime is far from trivia l 
matter (Bernardeau 1995; Gaztahaga & Bernardeau 1997). 



Smaller scales could be penet rated as well by usin g either ex- 
tended perturbation theory (^ojomb^^j/^^97 ) , or second 
order perturbation theory ( Scoccimarro, & Frieman 199q; 



ically small bias formally b = 7/3£ ; , because the first term 
in the Taylor expansion does not describe the bias accurately 



scoccimarro, fc Frieman 1996^ in combination with expand- 
ing the bias correction further than leading order. Note that 
smoothing in this theory precedes biasing as in FG93. While 
these two operations do not commute, according to David 
Weinberg (1998, private communication) the effect of this 
appears to be negligible in simple models on large scales, 
although it could lead to problems for some exotic biasing 
schemes. Based on the group properties of the bias func- 
tion (FG94) it is possible to apply Equations 1-2, and their 
above mentioned possible generalizations, for biasing be- 
tween cl usters and galaxies, or between different types of 



when the variance is large. 



galaxies (Kaiser 198l|; [Bernardeau fc Schaeffer 1992 Szapud 



Szalay 1993b; Bernardeau 1996) to predict the enhance 



Besides systematic errors, which are expected to be 
quite small for the SDSS and 2dF, redshift distortions could 
pose a problem. For the two-poi nt function t he theory of 
redshift distortions is well known (Kaiser 1987), and can be 
corrected for. For Q3 the effects o f redshift distortions were 



calculated by ( Hivon et al. 1995 ), an d they found it to b e 
negligible on large scales; simulations ( Szapudi et al. 19981 ) 
confirm this on scales when the variance is small, which is 
anyway the limit of applicability of the theory described 
above. It is safe to assume that the same holds for Q21 as 
well, since it is a third order cumulant with slightly different 
smoothing than Q3. Thus redshift distortions do not consti- 
tute a real problem for applying the theory of the previous 
section. The light cone effect is even less of a problem than 
redshift distortions: higher order moments of the distribu- 
tion can be corrected for it, and in any case it is not expected 
to be important for the next generation of wide field surveys 
where the median redshift is z ~ 0.2 (Matsubara, Suto, & 



ment of the two-point correlation function on large scales 
from counts in cells. The mentioned group property enables 
more then one step as well, i.e. between mass to galaxies to 
clusters, etc. Note that all the previous considerations were 
based on local bias models. It is possible to construct more 
general models which are stochastic and/or non-local in na- 
ture. These possible improvements and generalizations, as 
well as a full sca le demonstration of the technique on mock 
galaxy catalogs (Cole el al. 1998) with different biases will 



Szapud ^ 1998| ). 

Practicality of measuring these key quantities was 
demonstrated by several authors: correlation functions, and 
cumulants w ere estimated in numerous gal a xy catalogs and 
simulat ions ( Peebles 198C |_ 3aztana^a^F992j^zapuai j _Szah^ 



& Bosc fian 1992t [Meiksin, Szapudi, fc Szalay 1992j |Bouchet 



et al. 1! 93 j Gaztahaga 1994 : Szapudi et al. 1995; Colombi et 



Szapudi. Meiksin, &i Nichol 1996| ), cumulant corre- 



al 1991 

lators were determined in the APM catalogj Szapudi fc Sza- 
lay 199_|7a|), and recently in simulations ( [Vlunshi fc Melott 
1998n71nfact, the slope —7 in the previous equations is pre- 
cisely the slope of the second moment of fluctuations in the 
same cells in which the cumulants and cumulant correlators 
are measured. Therefore in practice, only counts in cells and 
bicounts (or joint counts) in pairs of cells have to be esti- 
mated for obtaining all required quantities. Thus, there is no 



be presented elsewhere. 

In summary, a novel method has been presented to de- 
termine the bias parameter b, as well as the higher order 
bias parameters. It is simple to use and will become a prac- 
tical possibility in the near future when the new generation 
of wide field surveys, the SDSS and 2dF, come on line. The 
technique involves the estimation of (bi)counts in cells, and 
from them the slope of the correlation function —7, and two 
third order moments: Q3, and Q21. The bias parameter can 
be obtained explicitly as b — 7/6(Q2i — Q3) in the context of 
gravitational instability with Gaussian initial conditions. It 
was demonstrated that these quantities can be measured in 
the near future, and iV > 3 connected moments will provide 
additional constraints as internal consistency check. Thus a 
unique determination of the bias will be a realistic possibil- 
ity with higher order statistics. Since velocity measurements 
yield Q ' 6 /b, this method can be combined with bulk flow 
measurements to yield f2 as well. The general theory pre- 
sented here can be used to describe bias between different 
species of galaxies and clusters as well. 
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